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Abstract. We found 25 galaxies with probable tidal tails 
in the Hubble Deep Fields North and South at z— 0.5-1.5. 
General characteristics of the selected tidal features are 
very close to characteristics of tidal tails of local interact- 
ing galaxies. Using objects with 0=0.5-1.0, we found that 
volume density of galaxies with tidal tails changes with 
z as {1 + z)'^^^. Therefore, we estimated the rate of close 
encounters between the galaxies of comparable mass lead- 
ing to the formation of extended tidal structures. If this 
rate reflects the merger rate, our data support a steeply 
increasing merger rate at 2; '--^l and consistent with zero 
curvature universe. 
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1. Introduction 

The redshift dependence of the interaction and merger 
rate is an important test of the current models for the for- 
mation and evolution of galaxies. It is now well established 
that galaxy interactions play a major role in galaxy for- 
mation and evolution (e.g. Schweizer 1998, Combes 1999 
for recent reviews). Toomre (1977) demonstrated that the 
universe's higher density in the past (cx(l-l-z)'^ ) suggests 
a higher past merger rate, increasing back in time as 
i^'^/s ^.^yj^jj tinie t) if the binding energies of binary galax- 
ies had a flat distribution. If the galaxy merger rate is 
parametrized in the power-law form ex (1 -I- z)™, then the 
exponent has been found to be m = 2.5 from Toomre's 
(1977) approach (assuming il = 1). Statistics of close 
galaxy pairs from faint-galaxy redshift surveys (e.g. Yee 
& Ellingson 1995, Le Fevre et al. 1999) and morpholog- 
ical studies of distant galaxies support a large value of 
the exponent m for z < 1. For instance, Abraham (1998) 
concluded that current best estimates for the merger rate 
are consistent with rn « 3. Preliminary studies of distant 
peculiar objects representing distinct results of interac- 
tions/mergers (collisional ring galaxies, polar-ring galax- 
ies, mergers) also support m > 3 (Lavery et al. 1996, 
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Reshetnikov 1997, Remijan et al. 1998, Le Fevre et al. 
1999), although statistics are still insufficient. Many other 
surveys, including IRAS faint sources, or quasars, have 
also revealed a high power-law (e.g. Warren et al. 1994, 
Springel & White 1998). However, some recent works have 
suggested a moderate (m < 2) (e.g. Neuschaefer et al. 
1997, Wu & Keel 1998) or intermediate (m - 2 - 3) 
(Burkey et al. 1994, Im et al. 1999) density evolution of 
merging systems with z. 

From analytical formulation of merging histories (e.g. 
Carlberg 1990, 1991; Lacey & Cole 1993), it is possi- 
ble to relate the dark haloes merger rate to the parame- 
ters of the universe (average density O, cosmological con- 
stant A). The merging rates for visible galaxies should 
follow, although the link is presently not well known 
(Carlberg 1990, Toth & Ostriker 1992). Theoretical mod- 
els based on Press-Schechter formalism (Carlberg 1990, 
1991) predict a redshift evolution of the merger rate with 
m ex J7°'*^(l — A)"'^-^-'^ (the exponents must be some- 
what changed if the average halo mass decreases with z 
- Carlberg et al. 1994). This conclusion is confirmed by 
numerical simulations within the CDM scenario - m = 4.2 
(fi — 1) and m — 2.5 (fJ — 0.3) for z < 1 (Governato et 
al. 1997). 

Tidal tails originate in close encounters of disk galax- 
ies (e.g. Toomre & Toomre 1972). The purpose of this 
note is to show that statistics of galaxies with extended 
tidal tails (tidal bridges have, on average, fainter surface 
brightnesses - Schombert et al. 1990) is a useful tool to 
study evolution of interaction rate to z ^ 1. The simu- 
lations by Hibbard & Vacca (1997) - they showed that 
extended tidal features remain readily visible in the long 
exposures typical of the Hubble Deep Fields out to 2: ~ 1 - 
is the theoretical base for our work. We found that current 
statistics of such objects in the North and South Hubble 
Deep Fields (HDF-N and HDF-S correspondingly) leads 
to TO « 4. (Preliminary results based on the HDF-N only 
are presented in Reshetnikov 1999 ~ Paper I.) 
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Table 2. Mean characteristics of interacting galaxies 
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Fig. 1. The /-band images of the HDF-S galaxies with 
tidal tails. The size of the panel is 19" x 19" . 



2. Sample of galaxies 

We used the deepest currently available deep fields (HDF- 
N - Wilhams et al. 1996 and HDF-S - WiUiams et al. 1999) 
to search galaxies with extended tidal tails. From detailed 
examination of the fields in the F814W filter (hereinafter 
referred to as /), we selected more than 70 tailed objects. 
Careful analysis of their images in combination with the 
redshift data enabled us to distinguish 25 galaxies with 
z = 0.5-1.5 (12 objects in the HDF-N and 13 in the HDF- 
S). Galaxies with tidal tails in the HDF-N are described 
in detail in Paper I. Here we present the data for the 
HDF-S objects (Fig.l). (Our statistics of galaxies with 
tidal structures are in general agreement with van den 
Bergh et al. (1996) data on morphology of galaxies in the 
HDF-N. van den Bergh et al. classified 20 galaxies with 
21 < / < 25 as objects with probable tidal distortions in 
the HDF-N.) 

General characteristics of the galaxies are summarized 
in Table 1. The columns of the table are: galaxy iden- 
tification, / band apparent magnitude, photometric red- 
shift (there are no published spectroscopic redshifts for 
the sample galaxies) . All the data are taken from the web 
site of the HDF-S group at SUNY, Stony Brook (Chen et 
al. 1998). In the fourth column we present the absolute 
magnitude in the rest-frame B band calculated according 
to Lilly et al. (1995) as: 

MB=/-51og(DL/10pc)+2.51og(l+z) + (B - 4)+G.17, 
where 



Dl 



Hoq'o 



[qoz + {qo - l){V2qaZ + 1 - 1) 



Parameter 


HDF North 


HDF South 


Number of galaxies (all) 

z =0.5-1.0 

z =1.0-1.5 

< z > 

< Mb > 

< iiB >(tail) 

^tail 1 ^galaxy 


12 

7 

5 

0.94±0.30 

-18.0±1.3 

23.7±1.0 

1.4±0.7 


13 

7 

6 

0.93±0.28 

-18.9±1.5 

23.1±1.0 

1.2±0.5 



(1) 



- luminosity distance (A = 0), Hq - the Hubble constant 
(75 km/s/Mpc), go - deceleration parameter (go=0.05), 
B — Iz - /c-correction color (we used correction for Sbc 
galaxy), and term 0.17 translates AB magnitudes into 
standard B. 

In Table 2 we compare mean characteristics of the 
tailed galaxies in two deep fields. As one can see, both 
samples are consistent within quoted errors. 

3. Characteristics of tidal tails 

To be sure that our selected galaxies possess true tidal 
tails, we performed photometric measurements in the I 
passband using circular apertures centered on the bright- 
est regions of the suspected tails. For the measurements, 
we retrieved the HDF-S images (version 1) from the ST Scl 
web site and processed them in the ESO-MIDAS environ- 
ment. The results of these measurements are summarized 
in Tabic 1 (column 5). The observed surface brightness of 
the tails has been converted to a rest-frame B by applying 
the cosmological dimming term and a fc-correction color 
term: ^l{B) = n{I) - 2.b\og{l+zf + {B ~ h) + 0.17 
(Lilly et al. 1998). General photometric characteristics of 
the local tidal tails are close to those for late-type spi- 
ral galaxies (Sb-Sc) (Schombert et al. 1990, Reshetnikov 
1998) and we used color term for Sbc galaxy (Lilly et al. 
1995) in our calculations. The results are presented in the 
last column of Tabic 1 and in Fig. 2. 

The mean rest-frame surface brightness of the tidal 
structures in the joint (HDF-N plus HDF-S) sample is 

< fXB >(tail) = 23.4±1.1, in full agreement with our re- 
sults for the local sample of interacting galaxies (obtained 
by analogous manner) - < jib >(tail) ~ 23.8±0.8 (Reshet- 
nikov 1998). 

Fig. 2 (top) presents the observed distribution of the 

< n{I) > values of the suspected HDF tails (dashed line) 
in comparison with the distribution for local objects in 
the B passband. The bottom part of the figure shows con- 
verted to a rest-frame, B distribution for HDF tails. It is 
evident that tails of distant galaxies demonstrate a dis- 
tribution of < IJ,(B) > values close to that for local in- 
teracting galaxies. Fig. 2 illustrates clearly the influence of 
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Table 1. Characteristics of galaxies in the HDF-S 



Name 


/f814W 


Zph 


Mb 


^/(tail) 


/is (tail) 


SB-WF-0564-0515 


25.27 


1.01 


-17.4 


25.4 


23.1 


SB-WF-0954-0586 


22.63 


1.29 


-21.1 


24.9 


21.6 


SB- WF- 1360-0405 


24.26 


1.02 


-18.4 


25.1 


22.8 


SB-WF- 1836- 1844 


22.85 


0.91 


-19.4 


24.9 


23.0 


SB-WF-2505-1580 


24.04 


1.37 


-20.1 


25.4 


21.6 


SB-WF-2536-2934 


22.68 


1.00 


-20.0 


24.9 


22.6 


SB-WF-2580-1700 


22.62 


1.27 


-21.0 


26.2 


22.9 


SB-WF-2658-3739 


21.78 


0.47 


-18.7 


25.7 


25.1 


SB-WF-2782-3543 


23.42 


1.08 


-19.5 


24.7 


22.2 


SB-WF-2962-2059 


24.20 


0.51 


-16.5 


25.1 


24.4 


SB-WF-3477-3723 


25.79 


0.95 


-16.7 


25.8 


23.7 


SB-WF-3579-0685 


21.33 


0.58 


-19.7 


24.9 


24.0 


SB-WF-3816-0752 


24.54 


0.69 


-16.9 


24.4 


23.1 



observational selection on the recognition of tidal struc- 
tures - we are able to detect relatively faint tails among 
the galaxies with z = 0.5 — 1.0 but among z ~ 1.0 — 1.5 
objects we can see only very bright tails. Therefore, our 
sample of galaxies with extended tidal tails is sufficiently 
incomplete for z > 1. Thus, objects with z = 0.5— 1.0 will 
give a more reasonable estimation of m in comparison with 
the total sample. 

4. Density evolution 

The resemblance of morphological and photometric char- 
acteristics of suspected tidal tails of the HDFs galaxies 
with local objects allows us to use them to measure possi- 
ble change with z of volume density of galaxies with tails 
(and, therefore, the rate of close encounters leading to the 
formation of extended tails). 

The co-moving volume element in solid angle dfl and 
redshift interval dz is 



dV 



-(1 



^y 



^■dndz, 



E{z) 



(2) 



- photometric 



where D^ 

and E{z) = 

1993). The increase of the space density of galaxies with 

tidal tails we take in standard power-law form: 



distance (eq.(l)), 
1 for A = (e.g. Peebles 



n{z) = no (1 + zf' 



(3) 



where no — n{z = 0) - local volume density of such galax- 
ies. By integrating equations (2) and (3) we can find the 
expected number of objects within solid angle dO, and in 
required range of z. 

4-.1. Local density of galaxies with tidal tails 

We suppose that at the current epoch interactions and 
mergers accompanied by tail formation are almost entirely 
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Fig. 2. The distribution of surface brightnesses of tidal 
tails. Top: solid line - distribution for local galaxies in 
the B band according to Reshctnikov (1997), dashed line 
- distribution for objects in the HDF-N and HDF-S in 
the / filter. Bottom: distribution of the B band rest- frame 
surface brightnesses for HDF tails (calculated from ob- 
served / band distribution). Dashed line shows galaxies 
with z = 0.5 - 1.0. 

between bound pairs of galaxies (e.g. Toomre 1977). So 
we adopt that frequency of tidal tails among single ob- 
jects (mergers) and in groups, is significantly lower than 
in pairs. 

According to Karachentsev (1987), the relative fre- 
quency of galaxies with tails among the members of bi- 
nary systems is 94/974==0.10±0.01. The fraction of paired 
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galaxies in the local universe is not well determined. Var- 
ious strategies give results between 5% and 15%. For 
instance, local pairing fraction is 7%±1% according to 
Burkey et al. (1994), 6%-10% (Keel & van Soest 1992), 
14%±2% (Lawrence et al. 1989). The most intensive stud- 
ies lead to 12%±2% (Karachentsev 1987) and 10% (Xu 
& Sulentic 1991, Soares et al. 1995). Moreover, Xu & 
Sulentic (1991) found that the fraction of pairs is ap- 
proximately constant (10%) over the luminosity range 
-22 < Mb < -16 (see also Soares et al. 1995). Thus, 
we can adopt the value of 10%±5% as a reasonable esti- 
mate of the local fraction of binary galaxies. Therefore, 
the relative fraction of galaxies with tidal tails at z = is 
0.1x0.1=0.01±0.005. 

To find total density of galaxies in the nearby part 
of the universe (z < 0.05), we considered the galaxy lu- 
minosity function (LF) according to Marzke et al. (1998). 
The adopted Schechter function parameters of the LF are: 
M* =-20.05, (f)*=5AxlO-^ Mpc^^ and a=-1.12 {Hq =75 
km/s/Mpc). By integrating LF from Mb = —15.4 to - 
21.1 (the range of absolute luminosities of galaxies with 
tails in the HDF-N and HDF-S), we found that total vol- 
ume density of galaxies is equal to 0.026 Mpc~^. Thus, 
no =0.01x0.026=(2.6±1.3)xl0-4Mpc-3. 

The total angular area within which we searched tailed 
galaxies in two HDFs is 10.4 arcmin^ or 8.8xl0~^ sr. 

4-. 2. Exponent m from tidal structures 

Varying exponent m, we can estimate the expected num- 
ber of galaxies with tidal features in the HDFs. In Fig. 3 
we present the results of calculations for two redshift in- 
tervals: 0.5-1.5 (total sample) and 0.5-1.0 (adopted cos- 
mology is A = 0, (7o = 0.05 or f7 = 2qa = 0.1 and Ho^75 
km/s/Mpc). As one can see, the total sample (25 objects) 
leads to m = 2.6. But this value must be considered as 
a low limit only due to strong underestimation of tidal 
tails at z > 1 (sect. 3). For the galaxies with z = 0.5 — 1.0 
(A'^=14) we obtain m = 4.0. Assuming Poisson error of 
N (\/]V=3.7), we have m = 4.0l!];5. Adding 50% uncer- 
tainty in the local space density no, we have obtained a 
final estimation of ?ti as 4.0to'9- (^^^^ ^^ note also that two 
potential sources of errors - underestimation of no value 
and omission of tailed galaxies in the HDFs - bias value 
of m in opposite directions and partially compensate each 
other.) 

The value of m depends on the adopted cosmologi- 
cal model. For A = 0, O = 1 we have m = 4.9 (z=0.5- 
1.0). In our calculations for the model with a cosmologi- 
cal constant and zero spatial curvature (rim=0.3, riA=0.7, 
Vt = rim + r2A=l) we used the analytical approximation of 
the luminosity distance Dl according to Pen (1999). In the 
framework of that model we have obtained m — 3.6lg'g. 

To obtain more realistic error estimation, we must take 
into account the possible luminosity evolution of galax- 
ies with redshift. Unfortunately, luminosity and surface 
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Fig. 3. The dependence of expected number of galaxies 
with tidal tails in two HDFs on exponent m for z=0.5-1.5 
(dotted line) and z=0. 5-1.0 (solid line). Horisontal lines 
show observed quantities. 



brightness evolution of peculiar and interacting galaxies 
is poorly constrained at present (e.g. Roche et al. 1998). 
Moreover, Simard et al. (1999) claim that an apparent sys- 
tematic increase in disk mean surface brightness to z ^ 1 
for bright {Mb < — 19) spiral galaxies is due to selection 
effects. Nevertheless, assuming that interacting galaxies 
undergo luminosity evolution A Mb = 1™ between z = 
and 1, we estimated that the value of m must be decreased 
by A m w 0.5: m = 3.5 for VL = 0.1, A = and m = 4.4 
for n = l, A = 0. 

5. Discussion and conclusions 

On the basis of analysis of the HDF-N and HDF-S im- 
ages we selected 25 galaxies with probable tidal tails with 
z=0.5-1.5. Integral photometric characteristics of the sus- 
pected tails are close to that for local interacting galaxies. 
Considering the subsamplc of tailed galaxies with z=0.5- 
1.0 (14 objects), we estimated that co-moving volume den- 
sity of such galaxies changes approximately as (1+z)'*. 
(Hence the volume density of tailed galaxies at z = 1 
is n(z = 1) = 4 X 10^^ Mpc^^ for q^ = 0.05.) Inclusion 
in the sample of the galaxies with tidal bridges does not 
noticeably change the value of the exponent (Paper I). 
Therefore, we estimated the change of the rate of close 
encounters leading to the formation of extended tails. If 
this rate reflects the merger rate, we have obtained evi- 
dence of a steeply increasing merger rate at z ^ 1. (Our 
result is related to field galaxies. The evolution in clusters 
might even be stronger than in the field. For instance, 
van Dokkum et al. 1999 found m = 6 ± 2 for the merger 
fraction in rich clusters of galaxies.) 

How does our estimation of m agree with values ob- 
tained by other methods? The recent surveys of the evo- 
lution of galaxy pairs with z are consistent with to w 3 
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Fig. 4. The dependence of the relative fraction of irregular 
and peculiar galaxies on redshift. The data for z=0.4 and 
0.8 are from Brinchmann et al. (1998); for z=0.54 are from 
Roche et al. (1998); the local value (z=0) is estimated from 
Marzke et al. (1998) LF. Dashed line - (l+z)^ relation. 



(see references in Abraham 1998). Evolution of the rate 
of interactions according to our data is characterized by 
close (within quoted errors) value of to. Direct analysis 
of the morphology of distant galaxies at z < 1 suggests 
a significant increase of the fraction of irregular and pe- 
culiar systems with redshift (Fig. 4). If interactions and 
mergers are responsible for the observed asymmetries of 
galaxies (e.g. Conselice & Bershady 1998), this increase 
can reflect the increase of the interaction rate with z. As 
one can see in Fig. 4, relative fraction of Irr/Pec galaxies 
changes in accordance with m — 4. Nairn et al.'s (1997) 
result (35%±15% of peculiar galaxies down to / = 24.0) 
agrees with to = 4 also. Many other observational surveys 
and numerical works indicate a large (>3) exponent m 
(Sect.l). Comparison with predictions of analytical and 
numerical works shows that current observational esti- 
mates of the merger rate favor a zero curvature (17 = 1) 
universe (e.g. Carlberg 1991, Governato et al. 1997). 

Our results indicate that further detailed statistics of 
galaxies with tidal structures will be a powerful tool to 
quantify the interaction and merging rates evolution. 
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